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ABSTRACT 

We study the dust properties of galaxies in the redshift range 0.1 < z < 2.8 observed 
by the Herschel Space Observatory in the held of the Great Observatories Origins Deep 
Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity 
(Lm) and dust temperature (Td us t) of galaxies are derived from the spectral energy 
distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS 
and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR 
luminosities and dust temperatures of local galaxies at z < 0.1 using AKARI and IRAS 
data in the field of the Sloan Digital Sky Survey. We compare the Ltr, — Td us t relation 
between the two samples and find that: the median Td us t of iiersc/id-selected galaxies 
at z >0.5 with Lir >5x10 10 L q , appears to be 2-5 K colder than that of AKARI- 
selected local galaxies with similar luminosities; and the dispersion in Tdust for high-z 
galaxies increases with Lir due to the existence of cold galaxies that are not seen 
among local galaxies. We show that this large dispersion of the Lir — Td us t relation can 
bridge the gap between local star-forming galaxies and high-z submillimctcr galaxies 
(SMGs). We also find that three SMGs with very low T dus t (< 20 K) covered in this 
study have close neighbouring sources with similar 24- /im brightness, which could lead 
to an overestimation of FIR/(sub)millimeter fluxes of the SMGs. 

Key words: galaxies: evolution - galaxies: formation - galaxies: general - galaxies: 
high-rcdshift - galaxies: starburst - infrared: galaxies 



1 INTRODUCTION 

Understanding star formation (SF) mechanisms in galax- 
ies and how star formation rate (SFR) of galaxies evolves 
through cosmic time are key issues in the study of galaxy 
formation and evolution. The SFR is closely related to dust 
properties such as temperature (Td us t), mass, opacity, emis- 
sivity, and spatial extent. In order to investigate dust prop- 
erties of galaxies, it is important to obtain data covering 
the complete infrared (IR) wavelength range that includes 
both the "Wien" and "Rayleigh- Jeans" sides of the peak of 
the IR spectral energy distribution (SE D). Thanks to the ad - 
vent of the Herschel Space Observatory (jPilbratt et al.ll2010l ) 
with its very wide wavelength coverage (70 — 500/uu), we are 
now able to study complete IR SEDs of high-z galaxies at 
z ~ 1- 3. 

The "Rayleigh- Jeans" side of the IR SED of a galaxy 
is crucial for the study of dust properties, but up to now 
available photometric data for high-z galaxies have been 
limited to a small number of wavelength windows and to 
small regions in the sky. Despite the modest wavelength 
coverage of IR SEDs for high-z galaxies, it was suggested 
that luminous infrared galaxies (LIRGs; 10 L© < Lm < 
10 12 Lq) and ultraluminous infrared galaxies (ULIRGs; 
1O 12 L < im < 1O 13 L ) at high redshifts may be colder 



than t heir l o cal counterpa r ts (e. g., Rowan-Robinson et al.l 



2004 120051: ISaiina et all l200fj; ISvmeonidis et al l l2009t 



Seymour et alj" 20ld: seealso Muzzin et all 201oh . For ex- 



ample, Svmeonidis et all |200£ ) found that the IR SEDs of 
high-z galaxies, on average, are peaked at longer wavelengths 
than those of local galaxies by comparing Spitzer 70-^im- 
selected galaxies at 0.1 ^ z < 2 with IRAS 60-/im-selected 
galaxies at z < 0.1. They also found that the peak of IR 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA, 
f E-mail: hoseong.hwangScea.fr 



SEDs for local galaxies is shifted to the shorter wavelengths 
as IR luminosity (Lir) increases, while the peak of IR SEDs 
for high-z galaxies is located at a wider wavelength range 
comp ared to local ga l axies and varies little with Lir. How- 
ever. iMagnelli etakl (|2009l ). using the Spitzer data in the 
fields of Great Observ atories Origins Deep Survey (GOODS; 
iDickinson et all 120031 ') and Far Infrared Deep Extragalactic 
Legacy survey (FIDEL), suggested that IR SEDs of high- 
z galaxies at 0.4 < z < 1.3 are not significantly different 
from those of local galaxies. Several studies have also been 
carried out which suggest that the appa rent change with 
redshift is mainly a selection effect (e.g.. [Pope et al]|2006l : 
IChapin et ai1l2010h . 

It certainly seems that at high redshifts there are 
much colder populations of galaxies than local galax- 
ies with similar IR luminosities [e.g., (sub) mill i meter 



galaxies (SMGs) , iBlain et al. 20021: Chapman et all 
' Pope et all 120061: iKovacs et all 120061: iHuvnh et all 



2005; 



2007 



Coppin et all 120081 : IClements et al.l 120081 ] The relation 



between IR luminosity and dust temperature has been 
used as a useful tool in order to understand the con- 
nection bet ween galaxy populations such as (U)LIRGs 
and SMGs jDunne et al.ll200d: IChapman et alll2003l. [2005I : 



Kovacs et all l200d: IChanial et all 120071; lYang et all J200 
Younger et"all 120091: iDve et all 120091: IChapin et all 120091 . 
20ld : IClements et a.1.1 120101 : iMagdis et al.ll2010al ). Interest- 
ingly, in the Lir — Td ust plane, SMGs form a separate locus 
from the (U)LIRGs, which can suggest a different origin be- 
tween the two. 

In this paper, we investigate the IR SEDs and the dust 
properties of high-z galaxies by taking advantage of the 
wide wavelength covera ge (70 — 500^m) of the Photodetector 
Array Camera (PACS: [Poglitsch et all l2010h and Spectral 
and Photometric Imaging Receiver (SPIRE; iGriffin et al.l 
120101 ) instruments onboard Herschel. In order to study how 
dust properties of high-z galaxies are different from those 
of their local counterparts, we construct a sample of galax- 
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ies at z < 0.1 that wer e observed by the AKARI telescope 
l|Murakami et al.l l2007f ). The AKARI all-sky survey data 
contain flux density measurements up to 160 fj,m, which 
can probe the long-wavelength side of the peak of IR SEDs 
of local galaxies in a way similar to that of Herschel for 
high-z galaxies. Throughout, we adopt h — 0.7 and a flat 
ACDM cosmology with density parameters Qa = 0.73 and 
Sl m = 0.27. 

2 OBSERVATIONS AND DATA 

GOODS-North (hereafter GOODS-N) was observed by Her- 
schel as part of the Guaranteed Time Key Programmes 
PACS Extragalactic Probe (PEPlJ) and Herschel Multi- 
tiered Extragalactic Survey 

(HerMEsQ; Oliver et al. 2010, 
in prep.). Source extraction on these PACS and SPIRE im- 
ages was performed at the prior positions of Spi tzer 24- /xm- 
selected sources, and d etails are described in iBerta et all 
( |2010h . iMagnelli etHl (|2010T > and iRoseboom et all (|2Qlof )7 
PACS measurements are above the confusion limit, and we 
used flux densities in PACS bands down t o 3cr limits of 3 an d 
5.7 mjy at 100 and 160 /im, respectively l|Berta et al.| [2010). 
SPIRE measurements are used down to 5a limits of 4.4, 4.8 
and 7.6 mjy at 250, 350 and 500 /im, respectively. It is noted 
that these measurements lie below the 1-cr SP IRE confusion 
limit of 5.8, 6.3, 6.8 mjy (|Nguven et alJl201Ch . However, this 
limit is a spatially averaged statistical limit which consid- 
ers that galaxies are homogeneously distributed in the field 
and all affected in the same way by close neighbours. In 
this study, we flag galaxies more "isolated" than others for 
which SPIRE flux densities can potentially be more robust 
by using the higher spatial resolution 24 fj,m images (to be 
explained in i|3}. 

By combining these catalogues with the existing multi- 
wavelength data, we made a band-merged catalogue of IR 
sources having flux density measurements at Spitzer MIPS 
24 and 70 fj,m, Herschel PACS 100 an d 160 (an, and H er- 
schel SPIRE 250, 350 and 500 /an (see lElbaz etaI]|201Cl for 
details). For 493 galaxies detected in at least one out of the 
five PACS/SPIRE bands, we computed the IR luminosity 
using the SED models of ICharv fc Elbazl l|200ll . CE01) by 
allowing renormalization of the templates. The SED fit was 
applied to the flux densities at A rcst 30/xm, and Lir for 
all 493 galaxies was calculated. 

In order to determine the dust temperatures, we fit the 
observational data with a modified black body ( MBB) model 
by fixing the emissivity p arameter to j3 = 1.5 (|Hildebrandl 
1 19831 ; iGordon etal]|201oh . We used only the galaxies that 
satisfy the following conditions: 

(i) there should be at least one flux measurement short- 
wards to the FIR peak (i.e. 0.55A P < A rcst < A p where A p 
is the rest-frame, peak wavelength of the best-fit CE01 SED 
model and A rcst is the rest-frame wavelength of the observed 
data) ; 

(ii) there should be at least one flux measurement long- 
wards to the FIR peak (i.e. A p ^ A rost ); and 



1 http:/ /www.mpe.mpg.de/ir/Research/PEP 

2 http: / /hermes. Sussex. ac.uk 



(iii) the FIR SED should be physical (convex, not con- 
cave), i.e. we reject galaxies with (Sioo or Sieo) ^ 5*250 
and S250 < (S350 or S500), or (S100 or Sieo) > S350 and 
S350 ^ 5500. 

The above criteria ensure a well-sampled SED around 
the peak of the FIR emission, and the wavelength cut of 
0.55A P is introduced to reduce the contribution of warm 
dust (i.e. emission from very small grains). For these galax- 
ies, we fit the flux densities at A rost ^ 0.55A P using the 
MBB model. We then select galaxies having flux measure- 
ments on both sides of the peak for the best-fit MBB model, 
and secure a final sample of 140 galaxies. Uncertainties of 
Lir and Td ust were computed by randomly selecting flux 
densities at each band within the associated error distri- 
bution (assumed to be Gaussian) and then re-fitting. It is 
noted that we checked how Td ust measurements of galaxies 
are affected if we remove some data points for the fit. We 
found that Tdust becomes systematically lower/higher if you 
use only data points shortwards/longwards to the FIR peak. 
However, Td us t does not change systematically even if you 
remove some data points at A rcst ^ 0.55A P for the fit as 
long as you have measurements on both sides of the peak 
wavelength. 

For local galaxies, we construct an IR catalogue by 
cross-correlating our IR sources with the galaxies in the red- 
shift catalogue as foll ows. We use the IR AS Faint Sources 
Catalog - Version 2 (jMoshir et al.l H"992l ) . which contains 
173,044 sources with flux density measurements at 12, 25, 
60 and 100 /im. We also use th e j4L".Aii!L/Far-Infrared Sur- 
veyor fFIS: lKawada" et al. 120071 ) all-sky survey Bright Source 
Catalogue (BSCQ) ver. 1.0 that contains 427,071 sources 
over the whole sky, with measured flux densities at 65, 90, 
140 and 160/im. Among the measurements at IRAS and 
AKARI bands, we consider only the reliabkjfl flux densities. 
For the redshift catalogue, we use a spectroscopic sample 
of galaxies in the Sloan Digital Sky Survey Data Release 
7 (|Abazaiian et all 120091 SDSS DR7) complemented by a 
photometric sample of galax ies whose redshift information 
is available in the literature l|Hwang et alj|2010h . 

By adopting the same method applied to high-z galax- 
ies, we compute Lir for the local galaxies by fitting the 
CE01 templates to the flux densities at A rost ^ 30/xm. For 
Tdust estimates, in addition to the conditions (i) and (ii) used 
for high- z galaxies, we use only sources detected at 140 or 
160 /mi. This condition ensures a fair comparison with the 
high-z samples as in both samples have flux density mea- 
surements longwards to the FIR peak. Since IRAS 60 /xm 
and AKARI 65 /xm (or IRAS 100 /xm and AKARI 90 /xm) 
partially overlap, we use only IRAS data (Seo and 5*100) for 
the fit to avoid over- weighting when both IRAS and AKARI 
flux densities are measured. Finally, we use only the galaxies 
having flux density measurements on both sides of the peak 
of the best-fit MBB model. 

Fig. [1] represents example SEDs for a high-z galaxy in 
GOODS-N and a local galaxy in SDSS with the best-fits 
CE01 and MBB models. In Fig.d we plot Lir and T dus t for 



3 http:/ /www. ir.isas.jaxa.jp/AKARI/Observation/PSC/ 
Public/RN/AKARI-FIS-BSC-Vl-RN.pdf 

4 Flux quality flags are either 'high' or 'moderate' for IRAS 
sources and 'high' for AKARI sources. 
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Figure 1. Example SEDs for (a) a high-z galaxy in GOODS-N 
with L m = 6.7 X 10 11 (L Q ), T dust = 37 K, and z = 0.64; (b) a 
local galaxy in SDSS with L m = 4.9 X lO n (L ), T dust = 35 K, 
and z = 0.02. The best-fit CE01 and MBB models are indicated 
by solid and dashed lines, respectively. 
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Figure 2. (a) L IR and (b) T dust for 140 galaxies in GOODS- 
N (circles) and 190 galaxies in SDSS (diamonds) vs. redshift. 
Galaxies with spectroscopic and photometric redshifts are filled 
and open circles, respectively. Galaxies hosting AGN are indicated 
by red symbols. 



140 high-z and 190 local galaxies as a function of redshift, 
which shows that the two samples are distributed over simi- 
lar ranges of Lir and Td us t. The high-z galaxies in GOODS- 
N are found at 0.07 < z < 2.74 with 2.3 x 1O 9 (L ) < L m < 
5.3 x 1O 12 (L ) and 19 (K)< T dust < 58 (K), while the local 
galaxies in SDSS are in the range 0.005 < z < 0.119 with 
2.7 x 1O 9 (L ) < Lir < 3.0 x 10 12 (L Q ) and 21(K)< T dust < 
49 (K). 



3 RESULTS 

We show the relation between Lir and T dust for local 
and high-z samples in Figs. [3] and [4] respectively. Fig. [3] 
shows the distribution of local galaxies in comparison with 
known ULIRGs and SMGs in the literature. We determine 
a smoothed median trend of T dust as a function of Lir for 
our local galaxies. Since the contribution of active galac- 
tic nuclei (AGN) to determining the median trend of T dus t 
for local and high-z galaxies could be different, we ex- 
clude them when determining the median. We call AGN 
those sources whose optical spectral types are found to be 
Seyferts, low-ionization nuclear emission-line regions (LIN- 
ERs) or composite galaxies in the emission li ne ratio dia- 
gram l|Baldwin et al.lll98ll ; iKewlev et al.ll2006l ). We classify 
only the galaxies at z > 0.04 as AGN due to the problem 
of the small (3") fix ed-size aperture for SDSS spectroscopy 
l|Kewlev et alJuOOd ). It is seen that T dust remains constant 
in the lower end of the luminosity range (Lir < 5 x 1O 1O L ) 
and increases as Lir increases. This trend has been in- 
dicated in the previous studies b y a shift of the peak in 
the IR SEP with increasing L m (ISoifer et all 1 19871. fl989l: 
Chapman et al. | |2003l; IChapin et al.ll2009l ; ISvmeonidis et all 



Lapr 

2009t "see also lAmblard et alJl201Ch . At L m > 10 12 (L 



our samples are smoothly connected to the locus of known 
ULIRGs at intermediate/high redshifts. We also note that 
changing AGN selection criteria has a small effect on our 
results. 

In Fig. H^a), we plot high-z galaxies along with the loci 
of ULIRGs, SMGs and local SDSS galaxies. The median 
trend of T dust for high-z galaxies is also determined by ex- 
cluding those galaxies hosting X-ray selected AGN as well as 
possibly blended sources. AGN have either Lx [0.5-8.0 keV] 
> 3xl0 42 ergs s _1 , a hardness ratio (ratio of the counts 
in the 2-8 keV to 0.5-2 keV passbands) greater than 0.8, 
Nh ^10 22 cm" 2 , or br oad/high- ionization AGN emission 
lines (Bauer et al. 2 0041 ). We regard a galaxy as 'blended' if 
it has at least one neighbouring source within ~one beam 
FWHM of the galaxy position at 24 jj,m and SPIRE images 
(20" for 24^m and 250 /xm, 27" for 350 ;im, and 46" for 500 
/im). For the 24- /xm image, the criterion is at least two neigh- 
bouring sources instead of one. In addition, we require that 
the flux density of a neighbouring source should be larger 
than 50% of that of the galaxy. Among 140 high-z galax- 
ies, we have 84 'cl ean' and 56 'blended' galaxies (see also 
iBrisbin et"aHl20ld for a purity index that is a measure of 
blending of a galaxy). Fig. |2£b) represents the distribution 
of T dus t for 'clean' and 'blended' galaxies including AGN. 
This shows that only 'blended' galaxies have several cold 
galaxies with T dust ~ 20 K, while 'clean' galaxies do not (to 
be discussed in detail in £14. 2p . 

In Fig. |4ja) , the median trend of T dust for high-z galax- 
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Figure 3. Td ust vs. Lir for galaxies in SDSS. Galaxies hosting AGN are indicated by yellow symbols. The thick dashed line is a smoothed 
median trend of Td ust for local SDSS galaxies by excluding those with AGN, and the dot-dashed li nes are its envelope that includes 
90% of the galaxies above and below the me dian. The known l ocal infrared gala xies (IRGs) (D00: [ Dunne et a.1.1 1200CT) . SMGs (C05: 
IChapman et alj|2005l. K06: lKovacs et al-llioO^) and ULIRGs (Y07: lYang et alj2007t Y09: lYounger et alj|2009l) are plotted with triangles, 
squares and star symbols, respectively Among the SMGs in common between this study and C05, those having no neighbouring sources 
(clean) are denoted by large filled circles, while those possibly contaminated by neighbouring sources (blended) are denoted by large 
open circles. 
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Figure 4. (a) Td us t vs. Lip. for galaxies in GOODS-N. Galaxies hosting AGN are indicated by yellow symbols. The thick dashed line is 
a smoothed median trend of T(j ust for local galaxies adopted in Fig. \3\ and the dot-dashed lines are its 90% envelope. The loci of known 
SMGs and ULIRGs in Fig. [3] are plotted as regions filled by cyan and coral colour, respectively Clean and blended galaxies are denoted 
by filled and open circles, respectively. The thick solid line is a median trend of T(j ust for galaxies in the GOODS-N field and the dotted 
lines are its 90% envelope (filled with orange colour). Inclined, grey long-dashed lines indicate the AKARI selection function at the 
maximum redshift of local galaxies [z = 0.119). (b) Distribution of Td ust for high-2 galaxies. Clean and blended galaxies including AGN 
are denoted by hatched histograms with orientation of 45° (/ / with red colour) and of 315° (\\ with blue colour) relative to horizontal, 
respectively. 
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Figure 5. Deviation of Td ust (ATd us t) for high-z LIRGs from the 
median trend of local LIRGs vs. redshift. Filled and open sym- 
bols arc 'clean' and 'blended' galaxies, respectively. Star symbol 
represents the median value of ATd ust for 'clean' galaxies in each 
redshift bin. 



ies also appears fairly constant at low luminosities, but in- 
creases as I/ir. increases at higher luminosities. By compar- 
ing with local galaxies, we see that the median trend of 
Tdust for ifersc/ieZ-selected high-2 galaxies at z > 0.5 be- 
comes smaller than that for AKARI-selected low- 2 galaxies 
by 2 - 5 K at high luminosities (L IR > 5 x 10 10 L Q ). Note 
that this does not necessarily mean that high-2 galaxies are 
systematically colder than local galaxies with similar lumi- 
nosities (to be discussed in ij4.ip . It is also seen that Tdust 
scatter increases by about 5 K in the sense that the upper en- 
velope for high-z galaxies is similar to that of local galaxies, 
while the lower envelope for high-z galaxies stretches below. 
These results are consistent w ith the results based on the 
peak position of the IR SEDs (|Svmeonidis et alj 12009). In- 
terestingly, the high-2 galaxies with low Tdust fill the gap 
between local star-forming galaxies and high-2 SMGs in the 
Lxr — Tdust plane. It is also interesting that Tdust for AGN 
follows a similar trend as for star-forming galaxies, which 
indicates that the star formation of galaxies having AGN 
does not significantly differ from that of normal star- forming 
galaxies seen in the F IR regime, where SF dominates the 
IR emission of AGN dElbaz et all l20ld ; IShao et al l l201fj 
lHatziminaoglou et al.ll2010T l. 

To investigate how Td us t for high-2 galaxies evolves with 
redshift compared to local galaxies with similar Lir, in Fig. 
[5] we plot the deviation of Td us t for high-2 LIRGs from the 
median trend ((Tdust, local)) thick dashed line in Fig. [4^,) of 
Tiust for local LIRGs (AT dusl = T d u S t - (T d u S t,iocai)) as a 
function of redshift. This shows that ATd us t is negative, 
on average, which indicates that //ersc/iei-selected high-2 
LIRGs at 0.3 < 2 < 1.4 appear to be colder than AKARI- 
selected local LIRGs. ATdust changes little with redshift, but 
it is difficult to draw any strong conclusions due the small 
sample size. Since the deeper observation of GOODS fields 
will be conducted in PEP and in the GOODS-Herschel key 
programme (PI: D. Elbaz), the future Herschel data will help 
us to address this issue with a larger number of galaxies up 
to higher redshift. 

To investigate how the SEDs of high-2 galaxies are dif- 
ferent from the local SED templates of CE01, we plot the 
average SEDs (from UV to FIR) of cold and normal LIRGs 
at high redshifts in comparison with the local SED tem- 
plates of CE01 in Fig. [6] The high-2 galaxies are classified 
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Figure 6. Average SEDs of LIRGs in GOODS-N with spectro- 
scopic redshifts normalized by the luminosity at 80 fim. Observed 
data for 10 cold and 26 normal galaxies are plotted by filled and 
open circles, respectively. The average SEDs of the data points 
are plotted as solid lines. Colour-coded dashed lines are CE01 
templates normalized by the luminosity at 80 /im, and Lir, of the 
templates with the maximum and the minimum Lj R are labelled. 



into three groups by comparing their Td ust values with the 
median trend of Td us t for local galaxies (thick dashed line in 
Fig. : cold (those with Td us t values less than the lower 
envelope of local galaxies); normal (those with Td ust val- 
ues within the envelopes of local galaxies); and warm (those 
with Tdust values larger than the upper envelope of local 
galaxies). Fig. [6] shows that the peak wavelength of the FIR 
SED for cold galaxies is larger than that for normal galax- 
ies as expected. The average FIR SED of normal galaxies 
appears to match well with the CE01 template of similar lu- 
minosity, but that of cold galaxies is matched with roughly 
80 times less luminous template, even if the median Tir of 
normal and cold galaxies are similar ((Lir/Lq) ~ 5 x 10 
and 4 x 10 1 L© for normal and cold galaxies, respectively). 

We wish to emphasize that the longer-wavelength side 
of the peak of the FIR SED for cold galaxies could be af- 
fected by blending at long wavelengths. Since the beam 
FWHM gets larger as wavelength increases, the flux den- 
sities at longer wavelength can be systematically overesti- 
mated if there is more than one source within the resolu- 
tion limit. In addition, galaxies having neighbouring sources 
seem to have low Td us t compared to 'clean' galaxies, as seen 
in Fig.Ufb), which suggests that we should be cautious when 
we measure Tdust using multi-wavelength data with differ- 
ent resolutions. We have calculated the fraction of blended 
galaxies among normal and cold galaxies, and find values of 
34±6% and 58±9%, respectively. We have tried to eliminate 
such problem in our samples by identifying blended galax- 
ies, and by using the FIR source catalogue based o n prior 
extraction ijBertaet alj|20ld : iRoseboom et alj|2010h . How- 
ever, since we may not have completely deblend all sources 
and do not account for galaxies whose neighbours are less 
bright than 50% of the target galaxy, there could be a re- 
maining contribution to the FIR flux of cold galaxies due to 
this effect. 

Interestingly, the cold galaxies are brighter than the 
normal galaxies in the optical bands. This could imply that 
there is a correlation between Tdust and the optical proper- 
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ties such as galaxy morphology, colour, or size. Thus, the 
study of the relation between optical properties and Td us t 
may provide us useful hints for the origin of the cold galax- 
ies (Le Floc'h et al. 2010, in prep.). 



4 DISCUSSION 
4.1 Selection Effects 

The AKARI data contain flux density measurements up to 
A = f60 fim (A rcst = 157 /im with a median redshift of 
z — 0.015 for local samples), while Herschel data have mea- 
surements up to A = 500 /im (A rcs t = 330 /im with a median 
redshift of z = 0.5 for high-z samples). Therefore, the Td us t 
difference between the two samples seen in Fig. |4ja) could 
be affected by the differ ence in selection effects. For exam- 
ple, [chiii^nD (2010) found results similar to this study, 
in the sense that the dust temperatures of L'LylS'T-selected 
(250—500 /im) high-z galaxies at z < 3 are systematically 
cooler than IRAS 60-/im-selected local galaxies. When they 
account the difference in selection effects between the two 
samples, they find no evidence for evolution of Td ust . On 
the other hand, if we consider ZZersc/ieZ-selected local galax- 
ies (z < 0.1) observed in the Herschel- ATLAS key pro- 
gramme (|Eales et al.ll2010T l. we can find cold galaxies with 
T duat 5, 20 K that are n ot seen in our local sample (see Fig. 2 
in lAmblard etafll201Ch . This seems to support the idea that 
the different selection effects between local and high-z galax- 
ies can explain the Td us t difference. However, note that there 
are few local galaxies with Lir ~ 10 11 — 1O 12 L0 (z < 0.1) in 
the sample of lAmblard et all (|2010l . see their Fig. 3), which 
are crucial for understanding the Td us t difference between 
local and high-z galaxies (see Fig. [4^i). Since a wider area 
will eventually be covered by the Herschel- ATLAS program, 
the upcoming data will help us to address this issue. 

To investigate the effect of the different selection ef- 
fect (i.e. observed wavelength and detection limit) on the 
difference of Td ust between local and high-z galaxies, and 
on the existence of cold galaxies (Td us t < 27 K and Lir > 
3 x 1O 11 L0) that are not seen in the local sample but seen 
in the high-z sample, we made a following experiment: we 
move the best-fit SEDs of high-z galaxies to local universe, 
and check whether they can be observed with AKARI de- 
tection limit used for local galaxies in this study. If we move 
all 140 high-z galaxies in Fig. |4]to z = 0.015 that is a me- 
dian redshift of local galaxies, the expected flux densities 
at 140 and 160 /im for more than 99% galaxies are above 
the minimum values of flux densities for our local samples 
(£140 = 1-9 and Sieo = 0.25 Jy), which means that no galax- 
ies would be missed with the AKARI detection limit. This 
confirms that the existence of cold galaxies at high redshifts 
that are not seen in the local sample, is not simply because 
of the different selection effect. 

For the extreme case, if we move high-z galaxies to 
z = 0.119 that is a maximum redshift of local galaxies, 67% 
of the galaxies would be observed. The effect of the AKARI 
detection limits in the Td ust — Lir plane is shown as a thick 
long-dashed line in Fig. [4] Galaxies lying on the right of this 
line would be detected. The steep slope of the line, indicates 
that the AKARI detection limit does not affect Td ust distri- 
bution significantly. Moreover, cold galaxies with Tdust < 27 



K and Lir > 3 x 10 1 L0 are free from this selection effect. 
This again confirms that the existence of cold galaxies at 
high redshifts is not because of the difference in selection 
effects. 

To check how Td us t measurements for our local galax- 
ies become different if we had flux measurements at longer 
wavelengths (A > 160/im) like for high-z galaxies, we have 
re-estimated T dust for 44 galaxies fo und in common between 
this study and iDunne et all Q2000, D00) by combining our 
flux densities with 850 /im observations from D00. We found 
that the two estimates agree very well with a median differ- 
ence of 0.2 K (rms scatter of 1.8 K). In addition, when we 
compare the observed 850 /im flux density in D00 with the 
expected flux density from our best-fit MBB model without 
850 /im data, two flux densities agree well within the errors. 
These also indicate that Tdust estimates with AKARI bands 
are not biased toward high Tdust- 

4.2 Cold Galaxies at high redshifts 

Tdust for galaxies is a function of total SFR per unit 
dust mass, the dust emissivity, and the geometry, which 
are closely related to the global star formation efficiency 
(SFE) . This SEE is known to be cont rolled by the gas den- 
sity (|Schmidtl [l959l ; iKennicuttl ll99Sft . which is connected 
to the spatial extent of star-forming regions. Therefore, 
the existence of cold galaxies at high redshifts may im- 
ply that the spatial distribution of dust in these galaxies 
is more extended than that in local galaxies with similar 
Lir. This is consiste nt with CO observational results for 
some high-z galax i es dTacconi et al.ll2006l: llono et al.l 



2009 



2010 



iDaddi et aLI l2010l: iTacconi et al.l 120101 : iGenzel et aL. 
see also iKaviani et al.ll2003l ). On the other hand, the exis- 
tence of cold galaxies at high redshifts may indicate large 
dust masses for these galaxies compared to local galaxies. 
If we assume similar gas-to-dust ratios for local and high-z 
galaxies at a given Lir, the detection of a large gas content 
in high-z g alaxies compared to l ocal galaxies could support 
our results (|Tacconi et al.ll2006l . l20irIlDaddi et alj|2010h . 

Interestingly, the cold galaxies at high redshifts (0.3 < 
z < 1.4) in Fig.|4ja), fill the gap between local star-forming 
galaxies and high-z SMGs, c onnecting the two populations 
(see also lMagdis et al]l2010bl ). It is noted that previous stud- 
ies on Tdust measurement of SMGs were based on only a few 
photometric data points that cover only long-wavelength 
side of the thermal SEDs. Therefore, it is necessary to 
re-estimate Td us t for SMGs using Herschel data to check 
how the locus of SM Gs in the Lir — Td ua t plane changes 
l|Magnelli et al.ll2010l ; Chanial et al. 2010, in prep.; Chapman 
et al. 2010, in prep.) . There are 13 SMGs in common be- 
tween this study and lChapman et al. I l|2005h . Among them, 
we have three SMGs with Td us t measured in this study (Td us t 
for the other ten galaxies are not measured because they did 
not fulfil the selection criteria in due to the low signal-to- 
noise ratios in some bands) , and show them as squares in Fig. 
|4ja) . They are found to be in the locus of known SMGs, indi- 
cating that our measurements are consistent with those with 
ground-based submillimeter studies. However, we can not 
determine T d ust for SMGs at 3x 1O 1O L < L m < 7x IO^Lq, 
which is important for understanding the connection be- 
tween local galaxies and SMGs (discussed in lMagnelli et all 
120101) . 
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For these 13 SMGs, we have checked the 'clean' flag 
defined in this study, and show them in Fig. [3] there are 
seven 'clean' and six 'blended' galaxies. The three SMGs 
with Tduat values measured in this study, are found to be 
'clean' galaxies. However, other three SMGs with Td ust < 20 
K in Fig. [3] are found to be 'blended' galaxies. We have 
checked whether or not the redshifts of neighbouring sources 
are different from these SMGs, and found that the redshifts 
of close neighbours for two SMGs (J123636.75 + 621156.1 
and J123651. 76 + 621221. 3) are not similar to their compan- 
ion SMGs. This indicates that the low Td us t of SMGs could 
be caused by an overestimation of the FIR/(sub)millimeter 
fluxes due to blending problems. One SMG (J123721.87 + 
621035.3 at z ~ 0.978) has a close neighbour at z ~ 0.969 
that is separated by 8.6" (~ 70 kpc), which suggests that 
galaxies in an interacting pair cou ld be colder than isolated 
galaxies (e.g.. iTacconi et al.ll2006r i. This issue - the relation 
between the galaxy merging stage and Td us t - needs to be 
investigated with a large number of galaxies. 
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